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ABSTRACT 



Aims. In this letter we characterise IRAS 12556-7731 as the first lithium-rich M-type giant. Based on its late spectral 
type and high lithium content, and because of its proximity in angular distance to the Chamaeleonll star-forming 
region, the star was misclassified as a young low-mass star in a previous work. 

Methods. Based on HARPS data, synthetic spectral modelling, and proper motions, we derive the astrophysical param- 
eters and kinematics of the star and discuss its evolutionary status. 

Results. This solar-mass red giant (T e g = 3460 ± 60 K and \ogg = 0.6 ± 0.2) is characterised by a relatively fast rota- 
tion (vsini ~ 8km s _1 ), slightly subsolar metallicity and a high-lithium abundance, A(Li)= 2.4 ± 0.2dex. We discuss 
IRAS 12556-7731 within the context of other known lithium-rich K-type giants. Because it is close to the tip of the red 
giant branch, IRAS 12556-7731 is the coolest lithium-rich giant known so far, and it is among the least massive and 
most luminous giants where enhancement of lithium has been detected. Among several possible explanations, we cannot 
preclude the possibility that the lithium enhancement and rapid rotation of the star were triggered by the engulfment 
of a brown dwarf or a planet. 

Key words. Stars: late-type, low-mass, fundamental parameters, abundances; Stars: individual: IRAS 12556-7731 



1. Introduction 

The source IRA S 12556-7731 was first identified by 
iPrusti et all ()1992f ) in an IRAS survey of pre-main sequence 
(PMS) stars in the Chamaeleonll (Chall) cloud. The 
source was associated with a relatively bright (R «12mag) 
off-cloud st ar. A later an a lysis o f the star's near-IR mag- 
nitudes led lLarson et al.l (|l998l) to conclude that its IR 
colours are more similar to those of giants than those of 
PMS stars. The source was then classified as a young stel- 
lar object ca ndidate using selecti on criteria based on Spitzer 
mid-IR data (I Alcala et al. II2008D. but it lacks the IR-excess 
' (jLarson et alj Il998t lAlcala et al. 1 12008|) typical of young 
stars. Since lithium is efficiently destroyed by convective 
mixing in the interior of low-mass stars when the temper- 
ature at the bottom of the convective layer reaches about 
, 2.5xl0 6 K, the presence of strong Lil A6707.8A absorp- 
tion represents an important criterion for identification of 
low-mass young stars. Thus, a late-type spectrum and the 
detection of a strong Li I A6707.8 A absorption in t wo VL T- 
FLAMES/Giraffe spectra, have led ISpezzi et~aTI (|2008ft to 
assign a young stellar object nature to IRAS 12556-7731. 
Surprisingly, this star turned out to be over-luminous by 
about two orders of magnitude on the HR dia gram with re 



spect to other young stars in the Cha II cloud (ISpezzi et al 
2008). This result and the conclusion by lLarson et al 
(| 19981) that IRAS 12556-7731 could be a giant star, trig- 
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gered a new high-resolution spectroscopic study. We charac- 
terise the star as a lithium-rich M-type giant. Several spec- 
troscopic studies of K-type giants in the past two decades 
have provided evidence of lithium enhancem ents in the at- 
mosphere of these stars (jKumar et a,l.lf20TlI and references 
therein). We study IRAS 12556-7731 in the context of other 
well known lithium-rich giants and discuss possible mecha- 
nisms for its lithium enhancement. 



2. HARPS spectroscopy 

With the aim of ascertaining the nature of IRAS 12556- 
7731, two high-resolution (R~T10,000) spectra (range from 
3800 A to 6900 A), obtained in February 2-3, 2011 with the 
HARPS (High Accuracy Radial velocity Planet Searcher) 
spectrograph at the ESO 3.6 m telescope in La Silla 
(Chile), were combined. The strong Lil absorption line 
at A6707.8A was the first feature to be immediately 
confirmed, for which we measured an equivalent width 
of 550(±10)mA. Likewise, the Ha line is in absorption. 
Overall, a number of strong molecular absorption bands, 
mainly of titanium oxide, can be identified, and they con- 
firm it as a very cool object. The interstellar (IS) Nai D 
(AA5890,5996 A) absorption components, distinguishable 
from the photospheric ones, are clearly detected (see Fig.[TJ). 
The narrower NaiD wings of IRAS 12556-7731 with re- 
spect to HD 173740 (see Fig. [I} definitely exclude it as a 
main sequence star. 

The mean radial velocity of IRAS 12556-7731, as drawn 
from the HARPS spectra, is i?V=67.68±0.02kms _1 , 
which is far beyond t he range of the Chamaeleon 
region (~15±2kms- 1 ; [C ovmo et all I1997D . Likewise, 
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Table 1. Physical parameters of IRAS 12556-7731. 
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Fig. 1. Portion of the HARPS spectrum of IRAS 12556-7731 in the Na I D region. The M5 III best template (HD 175865) 
is overplotted (dotted line). A dwarf M star with nearly the same temperature (HD 173740) is also overplotted with a 
dashed line. The interstellar absorption components (IS) are also indicated. 



the proper motion components of the star, fx a cos 5= 
— 13.1±5mas yr -1 , fj,g= 10.7±5mas yr" 1 ! a s retrieved from 
the PPMXL Catalog (jRoeser et al.l l201fl . are not com- 
patible with thos e of the Chamaeleon young stars (e.g., 
iFrink et al.lll998l) . Therefore, we conclude that the kine- 
matics of IRAS 12556-7731 is inconsistent with those of 
young stars in the Chamaeleon region and that the star 
is unrelated to the star-forming cloud. The detection of the 
interstellar Nai absorption components, for which we mea- 
sure an RV\s ~ 14 ± 2kms~ 1 consistent with the radial 
velocity of Chamaeleon, is an indication that such inter- 
stellar components can be produced by the Cha II cloud it- 
self and that IRAS 12556-7731 must be located at a much 
greater distance. Because it is unrelated to the star-forming 
region, IRAS 12556-7731 is most likely a field giant star, far 
behind the Cha II cloud. This confirms the suggestion by 
lLarson et all (| 19981) . 

3. Physical parameters and lithium abundance 

Using the combined (S/N«30) HARPS spectrum, we ob- 
tained first estimates of the physical parameters by a pro- 
cedure outlined in detail in Appendix [AJ This procedure 
provides us with the best effective temperature, gravity, 
and metallicity by comparing the combined HARPS spec- 
trum with template spectra of real stars with well known 
parameters via a \ 2 minimization criterion. The derived 
stellar parameters come from a weighted mean of the pa- 
rameters of the 50 reference stars (10 per each spectral re- 
gion) that independently match the target spectrum most 
closely. The results are T cS = 3460 ± 60 K, log g= 0.6 ± 0.2, 
[Fe/H]= -0.05 ± 0.10, and vsi nz= 6.2 ± 3.0 km s " 1 . F rom 
the log g vs. T e s relationship bvlHoudash elt et al.l (|2000l) for 
solar metallicity giants, we estimate log g(3500 K) ~ 0.4, in 
fair agreement with the value derived from the x 2 minimiza- 
tion. 

3.1. Lithium abundance 

As pointed out by iGratton fc D'Antonal ([1989;), a high 
value (^500 mA) of the A6707.8 A Lil line equivalent width 



implies that the resonance doublet should be strongly sat- 
urated and that deviations from LTE conditions may be 
important. This requires the use of other lithium lines. The 
spectral range covered by HARPS also allows us to inves- 
tigate the A6103.6A Lil line. 

The Li abundance A (Li) was derived by spectral syn- 
thesis of the Lil 6103.6 A line using the MOOG code 
(|Sneden et al.lll973l ). which assumes LTE conditions, along 
with the GAIA model atmospheres (Brott & Hauschildt 
2010, priv. comm.). The hyperfine lin e structure has b een 
considered following the guidelines by IWahlgrenl (l2005|) to 
compute the atomic parameters. The atomic and molec- 
ular line lists were taken from the VALEQ database. The 
physical parameters derived from the x 2 minimization ex- 
plained above were used as starting values for the synthesis. 
We let metallicity, vs'mi, A (Li), and microturbulence ve- 
locity £ vary, while keeping T e ff and logg fixed. The low 
temperature of the star makes the spectral synthesis rather 
difficult, because of the uncertainty on continuum determi- 
nation caused by several overlapping molecular lines. The 
continuum was determined by a spline fit made after vi- 
sually selecting regions free of absorption lines. We esti- 
mated that the uncertainty on the continuum determina- 
tion, which is the main source of error on A(Li), is less than 
5%. The results are shown in Fig. [2] From our best fit, we 
obtained values of A(Li)= 2.4 ± 0.2, [Fe/H]= -0.1 ± 0.2, 
and Dsini = 9 ± 2 km s _1 for a microturbulence velocity 
£ = 0.8 ± 0.2 km s -1 . The [Fe/H] and vsini values are in 
good agreement with those from the % 2 minimization, so 
average values were computed. The final adopted physical 
parameters are reported in Table [TJ 

An attempt at spectral synthes is of the lithium lin e at 
6707.8 A with MOOG (using the iReddv et~aTl 120021 line 
list implemented within the VALD database) resulted in 
large residuals (c.f. Fig. 2]), mainly owing to non-LTE con- 
ditions, such as the effects of over- ionization, and line asym- 
metry due to the convective motions, which mainly influ- 
ence the A6707.8A L il line formation in very cool stars 
(jCarlsson et al.lll994h . Since the EW of this line could be 
measured with high accuracy, we could estimate the non- 

1 http:/ /vald. astro. univie.ac.at/~vald/php/vald.php. 
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LTE Li a bundance throug h extrapolating the curves of 
growth bv lLind et al.l f|2009T ) to a temperature of ^3500 K, 
which yield A(Li) ~2.5. 

We thus conclude that a reliable value for the lithium 
abundance of IRAS 12556-7731 is A(Li)=2.4 ± 0.2. As a 
result, the low temperature and high lithium abundance 
of the star make IRAS 12556-7731 the coolest lithium-rich 
giant known to date. As shown below, our determinations 
of temperature and gravity place the star close to the tip 
of the red giant branch (RGB) on the HR diagram. 
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Fig. 2. Portion of HARPS spectrum in the interval around 
the lithium A6103.6 A absorption line. The top panel shows 
the synthetic spectra for three values of A(Li) at fixed 
metallicity, while the lower one shows the synthetic spectra 
for three values of metallicity at fixed A (Li). The best fit is 
for a sligthly subsolar metallicity and A(Li)=2.4. 



4. Discussion and conclusions 

The stellar parameters of IRAS 12556-7731 indicate that 
the star is a lithium- rich M- giant. To investigate its 
evolutionary statu s in more de t ail, w e used the tracks 
and isochrones by iGirardi et al.l (|2000l ) for solar metallic- 
ity stars. A comparison of the position of IRAS 12556- 
7731 in the log g vs. logT c ff diagram (Fig. [3]) showed 
that the physical parameters are consistent with an age 
~10Gyr. Reasonable values for the stellar luminosity and 
mass at that age, which are consistent with logT c ff, are 
log(L/L Q ) =2.99±0. 15 and M/M Q =1.0±0.2, respectively. 
The position of IRAS 12556-7731 on the HR diagram is 
shown in Fig. El along with sever al lithium-rich stars, 
as compiled bv iKumar et al.l (|201lh . Most of the known 
lithium-rich giants are concentrated within a luminosity 
range of 1.3< log(L/L Q ) <2, with only a minority having 
log(i/L©) >2.2. Among these luminous K-type objects, 



HD 39853 is the coolest one; nevertheless, IRAS 12556-7731 
is about 450 K cooler, placing it among the least massive 
and most luminous lithium-rich giants known so far. Being 
close to the tip of the RGB, one should then consider three 
possibilities for the evolutionary stage of the star: i) giant 
branch ascent, i.e. H-shell burning phase; ii) AGB phase; 
or in) post He-core flash. Unfortunately, the CNO abun- 
dances and the 12 C/ 13 C carbon isotopic ratios, tracers of 
the degree of mixing that provide further constraints on the 
evolutionary status, cannot be derived because the spectral 
range of HARPS does not include the appropriate wave- 
length range. 
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Fig. 3. Upper panel. T he log<? versus logT e ff diagram. 
The IGirardi et al.l (|2000D isochrones for three different ages 
are overplotted. The black dot represents the position of 
IRAS 12556-7731. Lower panel. HR d iagram of several Li- 
rich giants. The evolutionary tracks bv lGirardi et al.l {2000) 
for several masses as labelled are represente d with contin- 
uous lines. The lithium-rich K-type giants (jKumar et al.l 
are over-plotted with open circles. The asterisk rep- 
resents HD 39853, while the position of IRAS 12556-7731 
is indicated by the big black dot. 



Assuming a distance cichaii = 178 pc, ISpezzi et al.l 
have derived a luminosity log(L/L Q ) ~1. By com- 



paring it with the luminosity estimated above, we conclude 
that IRAS 12556-7731 should be located at d «1.75kpc 
from the Sun. By adopting this distance and using the 
HARPS radial velocity and the proper motion compo- 
nents (cf. Sec. the resulting spatial velocity com- 
ponents of the star in a left-handed coordinate system 
are (U, V, W)=(+30, -121, +79) kms" 1 , where U, V, and 
W are directed towards the Galactic anti-centre, the 
Galactic rotation direction, and the North Galactic Pole, 
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respectively. According to the criteria of lOortl (fl9 26^PI. 
IRAS 12556-7731 can thus be considered as a high- velocity 
star. The Galactic latitude, distance, and kinematics of the 
star imply that it most likely belongs to the old thin-disk 
population, consistent with its almost solar metallicity. But 
what is the cause of its high lithium content? 

Three possible scenarions can be considered to ex- 
plain the high lithium abundance observed in IRAS 12556- 
7731: i) the star somehow preserved the original lithium in 
its atmosphere, ii) lithium has been regenerated in later 
evolutionary stages by t he Cameron-Fowler mechanisn0 
(|Cameron fc Fowlerll 19711 ) , Hi ) lithium was enhanced by the 
engulfment of a br own dwarf or planetary compa nion. 

As discussed in iFekel fc Balachandranl (|1993f) , the sce- 
nario of the preservation of initial lithium in red giants is 
very unlikely. They also note that high lithium abundance 
was not revealed in a sample of 200 F-type stars with masses 
between 1 and 2M Q , just evolved off the main sequence. 
As they point out, if anything, lithium preservation might 
eventually work in stars of the early-F and A types, that 
are more massive than IRAS 12556-7731. As a solar mass 
star, the initial lithium was most likely burned during the 
pre-main sequence phase. 

Depending on the stellar mass and luminosity, the pro- 
duction of lithium via the Cameron-Fowler mechanism su- 
possedly occurs both at the RGB lum inosity function bump 
([Charbonnel fc Balachandran [2000) and during the He- 
core flash ( Kumar et al.l 1201*1 ) . It is possible that lithium 
was synthesized in IRAS 12556-7731 by the Comeron- 
Fowler mechanism during the luminosity bump, but most 
likely it had already been destroyed as the star evolved up 
to the RGB. Although the lSackmann fc Boothrovdl (fl999) 
model predicts that a high lithium abundance may be pre- 
served all the way up the RGB, their parameterization 
requires very high mixing rates. It was assumed t h at ro - 
tation could provide such rates, but iPalacios et al.l {2006) 
find that a self-consistent model of rotational mixing cannot 
generate enough circulation to account fo r the mechanism 
work ing efficiently. It has been suggested (|de la Reza et al.l 
1996}) that the IR- excess observed in some lithium-rich gi- 
ants may be due to a dust shell possibly also generated 
via the Cameron- Fowler mechanism, b ut for IRAS 12556- 
7731 there is no ev idence of IR excess (jLarson et al.lll998t 
lAlcala et al. 1120081 ). On the other hand, it is unlikely that, 
assuming an AGB status for the star, lithium has been syn- 
thesized during the He-flash phase, because such a process 
shoul d work for stars wit hin a narrow mass range around 
2M (|Kumar et alJl201lh . A pure Cameron-Fowler mecha- 
nism should only provide 7 Li, with no 6 Li. An attempt at a 
simultaneous fit of the resonance and subordinate lithium 
line by including the 6 Li/ 7 Li isotopic ratio as a free param- 
eter improves the fit of the A6707.8 line (see Fig. [4]), but 
there is no way to get a good fit for values of 6 Li/ 7 Li less 
than 0.11. 

In the brown dwarf/planet engulfment scenario, the ac- 
creted matter would also result in a simultaneous enhance- 
ment of 6 Li, 7 Li, and Be. Unfortunately, our HARPS spec- 
trum does not achieve the wavelength range to investigate 
beryllium, but the suggestion that the 6 Li/ 7 Li isotopic ra- 
tio may be as high as 0.11 would support the accretion 



2 W + 10 ^SOkms" 1 and/or (U 2 + V 2 ) 1/2 >65km s" 1 . 

3 Conversion of 3 He to 7 Li by a-capture with 7 Be as a radioac- 



tive intermediary. 



scenario. Also, the large radius of IRAS 12556-7731 in com- 
parison with most of the known lithum-rich giants makes 
the planet engulfment scenario plausible, because planet in- 
gestion would be more likely to occur when a star evolves 
more in the RGB and achieves a larger radius. 

Finally, in discussing the case of IRAS 12556-7731, we 
have to consider that its projected rotational velocity, 
vsmi ~8km s _1 , is rather high in comparison with other 
lithium-rich giants. Several authors have argued that ac- 
cretion of a planet may expla in both lithium enhancement 
and rapid rotation in gian ts (jDenissenkov k, Herwid [20041 ; 
iCalsberg et al.l 120091 120101 and references therein) . Given 
its physical properties, we cannot rule out that such a pro- 
cess has led to the lithium enhancement and rapid rotation 
in IRAS 12556-7731. Some rapidly rotating giants can even 
become magnetically active and may be detected in X-rays 
([Guillout et all 120091 ), however IRAS 12556- 7731 was not 
detec ted in a ROSAT pointed observation ([Alcala et al. I 
120001) , 

With the observational data available so far, we cannot 
be conclusive about what process dominates the lithium 
enhancement in IRAS 12556-7731. Perhaps several mecha- 
nisms working at different times along the RGB phase of 
this star have contributed to its lithium enrichment. 
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Fig. 4. Comparison of the observed lithium A6707.8 A ab- 
sorption line with synthetic MOOG spectra for three values 
of A(Li). The 6 Li/ 7 Li isotopic ratio has been also included 
as a free parameter for the fit. The inset shows the com- 
parison of the results with and without inclusion of the 
6 Li/ 7 Li ratio as a free parameter (continuous and dotted 
lines, respectively). 




Fig. A.l. Portions of the HARPS spectrum of IRAS 12556- 
7731 around three TiO bands. The M5III best template 
(HD 175865) is overplottcd with dotted red lines. 



Appendix A: The ROTFIT procedure 

In order to derive the astrophysical parameters of IRAS 12556—7731, 
we used a \ 2 minimization procedure developed in the IDI0 environ- 
ment. Thi s procedure, which is a variation of the ROTFIT code de- 
scribed bv lFrasca et al.l fi 2QQ30 > provides us with the best match of the 
observed spectrum with a grid of high-resolution spec tra of real stars 
with well determined parameters (PASTEL ca talogue, ISoub iran et all 
2010) retrieved from the ELODIE archive jMoultaka et al.l l200lf h 
which is one of the largest datasets of high-resolution spectra, span- 
ning the wavelength range 4000-6800 A. Because the spectra are 
severely affected by molecular bands, it was difficult to perform a 
safe and homogeneous normalization to the local continuum for both 
IRAS 12556-7731 and the ELODIE templates. We thus preferred to 
analyse selected spectral regions, mainly around the strongest TiO 
bands (AA4950, 5166, 5450, 6158, 6650, 6680 A), normalizing the spec- 
trum with respect to the average stellar flux in a window of ~ 10 A 
blueward of the band- head (see Fig. lA.H . This allowed us to fully 
exploit the high sensitivit y to both tempe rature and gravity of the 
molecular bands (see, e.g.. lNeff et al.lll995t) . The NalD lines, which 
are very sensitive diagnostics of temperature and gravity, were selected 
as well (Fig.frj. At the end of this procedure, for each of the five spec- 
tral regions investigated, ten templates that best matched the HARPS 
spectrum of IRAS 12556-7731, based on \ 2 minimization, were cho- 
sen independently, and their known astrophysical parameters were 
used to compute weighted average values. Such average values were 
adopted as astrophysical parameters for IRAS 12556-7731. The results 
areT eft = 3460±60K, logg = 0.6±0.2dex, [Fe/H] = -0.05±0.10 dex, 
and vsin i = 6.2±3.0kms _1 . The errors include the la standard devi- 
ation on the average and the typical errors of the PASTEL astrophys- 
ical parameters added in quadrature. The latter are ±50 K, ±0.1 dex, 
±0.1dex, and ±0.5km s" 1 for T c ff, logg, [Fe/H], and vsini, respec- 
tively 



4 IDL (Interactive Data Language) is a registered trademark 
of IT Visual Information Solutions. 



